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jjj ; ABSTRACT 

With the advent of the COBE detection of fluctuations in the Cosmic Microwave Back- 
ground radiation, the study of inhomogeneous cosmology has entered a new phase. It is 
now possible to accurately normalize fluctuations on the largest observable scales, in the 
linear regime. In this paper we present a model-independent method of normalizing theo- 
ries to the full COBE data. This technique allows an extremely wide range of theories to 

t— i | be accurately normalized to COBE in a very simple and fast way. We give the best fitting 

^ ' normalization and relative peak likelihoods for a range of spectral shapes, and discuss the 
>- ■ normalization for several popular theories. Additionally we present both Bayesian and fre- 
quentist measures of the goodness of fit of a representative range of theories to the COBE 

q ! data. 

m ; 

Subject headings: cosmic microwave background - methods: statistical 

G\ '■ 

Introduction 

Classically, it has been standard practice to normalize models of large-scale structure 
at a scale of ~ 10/i -1 Mpc, using a quantity related to the clustering of galaxies (here 
the Hubble constant Hq = 100/ikms Mpc~ ) measured at the current epoch. Due to 
processing of the primordial spectrum and the large amplitude of the mass fluctuations 
.^J | which galaxies represent, this method of normalization requires assumptions about the 
^ ■ history of the equation of state for matter inside the horizon, the non-linear evolution of 
the density field and the processes of galaxy formation. One of the key uncertainties is 
the relationship between the observed structure and the underlying mass distribution in 
the universe. With the COBE DMR detection of Cosmic Microwave Background (CMB) 
anisotropies (Smoot et al. 1992), it has become possible to directly normalize the potential 
fluctuations at near-horizon scales, circumventing the problems with the 'conventional' 
normalization. 

In an earlier Letter we presented the normalization of the standard cold dark matter 
(CDM) model and a small range of variants (Bunn, Scott & White 1995). In this paper 
we extend this to a larger class of models, and present a means for normalizing a whole 
class of models to the COBE data in a computationally simple manner. Throughout we 
will use the 2-year COBE data (Bennett et al. 1994, Gorski et al. 1994, Wright et al. 1994, 
Banday et al. 1994, Gorski et al. 1995a) as released by NASA-GSFC. The normalization in 
this data differs from that in the data used by the COBE group by ~ 1/iK (K.M. Gorski, 
private communication; Gorski et al., 1995b). 
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As discussed in Bunn et al. (1995) and Banday et al. (1994) there is more information 
in the COBE data than just the RMS power measured. In other words, the COBE data 
cannot be reduced to a single number without a significant loss of information. One way 
to see this is to notice that there are ~ 90 eigenvalues of the signal-to-noise ratio with 
eigenvalue larger than 1. An alternative method is to note that the COBE data constrains 
the amplitude of the fluctuations over a range of scales, albeit a narrow range. 

In Fig. 1 we show the RMS power in a model which is fittetedd to the COBE data with 
one free parameter. The toy model we have chosen is the so called Sachs- Wolfe spectrum 
(Sachs & Wolfe 1967) which assumes that the observed temperature fluctuations come 
purely from the redshifts associated with climbing out of potentials on the last scattering 
surface (for a discussion see e.g. Peebles 1993; White, Scott & Silk 1994). Aside from the 
normalization of the model (which is fixed by COBE) there is one free parameter: the 
spectral slope, denoted n, with n = 1 corresponding to a scale invariant spectrum. We 
notice that the total power is not constant, showing that the normalization to COBE is 
sensitive to more than the total RMS fluctuations produced. Furthermore, the COBE data 
contain information on the "shape" of the power spectrum, which means that some theories 
are more likely than others, given the data. We will now introduce a simple method for 
using all of the information in the COBE data to normalize a wide class of models. 

Model independent analysis 

We have demonstrated that the COBE data contain more information than just the 
total power, which can be measured by (Q), (j(10°), cr(7°), band power, etc. However, it 
is still useful to have a method for normalizing a given model to the data, without having 
to use the full sky maps containing 6144 pixels at 3 frequencies in each of two channels. 
We present here a method which allows one to normalize a large class of theories (those 
which can be described by a power spectrum over the limited range of scales probed by 
COBE) to the data in a simple manner. 

To proceed we notice that theories with gaussian fluctuations (or fluctuations which 
are gaussian on COBE scales) can be specified entirely in terms of a power spectrum. 
In CMB studies this is usually expressed as the variance of the multipole moments, as a 
function of mode number, i.e. £{£ + l)Cg vs £, where the temperature on the sky has been 
expanded in spherical harmonics 

AT 

— (0,0) = $><m*WM) (1) 

£m 

and we define 

( a e m a e'm>) = Ci S££>5 rnm ' (2) 

For most theories the power spectrum is a smooth function. Following White (1994), 
we will write 

D(x) = £{£ + l)C t with x = log 10 £ . (3) 

We can perform a Taylor series expansion of D(x) about some fiducial point, which we 
shall take to be x = 1 (£ = 10). Many theories (see below) can be well approximated 
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by quadratic D(x) over the relevant range for COBE, roughly £ = 2 to 30, and so we 
present the normalizations and likelihoods of quadratic D(x). We choose to parameterize 
our quadratics by the (normalized) first and second derivatives at x = 1: D[ and D'{ where 

D(x) ~ £>i (l + Di(a; - 1) + ^(z - l) 2 ^) (4) 

(note that D[ is l/.Di times the derivative of D(x) at x = 1). The normalization is then 
given by quoting D 1; or C10, for each (D^, £)") pair, and the goodness of fit is quantified 
by the relative likelihood of that shape compared to a featureless, n = 1, Sachs- Wolfe 
spectrum. 

We compute the maximum-likelihood normalization for a grid of values of D[ and 
using the method described in Bunn & Sugiyama (1995) and Bunn et al. (1995). We 
combine the six publicly-available equatorial DMR sky maps pixel by pixel into a single 
map by performing a weighted average, with weights given by the inverse square of the 
noise level in each pixel. (We obtain negligibly different results when the two 31 GHz 
maps, which are more sensitive to Galactic contamination, are excluded. These maps have 
high noise levels, and are therefore automatically given low weights in the average.) We 
excise all pixels whose centers have Galactic latitude \b\ < 20° from the map, leaving 4038 
pixels. 

We then estimate the likelihood of getting this particular data set for each power 
spectrum. Rather than computing the likelihood directly, which would involve repeated 
inversions of a 4038 x 4038 matrix, we first perform a Karhunen-Loeve transform (Karhunen 
1947, Thierren 1992) to "compress" the data to a more manageable size. Specifically, we 
choose a set of basis functions fa defined on the portion of the sphere outside of the 
Galactic cut. (The manner in which these functions are chosen is described below.) We 
then compute the inner product of the data vector with each of these functions: 

3 

where Yj is the position of the jth pixel, AT(fj) is the corresponding temperature fluc- 
tuation in the data, and the sum runs over all pixels. If we assume that the temperature 
fluctuations are gaussian, then the projections Xi will be gaussian as well. We can therefore 
compute their likelihood in the usual way: 

L oc (det V)~ 1/2 exp (-^XiV^Xj) , (6) 

where Vij = (xiXj) is the covariance matrix, which contains contributions from the cosmic 
signal and the noise: 

V = V sig + Koise = (FY) (BCB) (FY) T + FNF T , (7) 

where = /*(?_,-), Y ifl = Y im (ri), C^ v = Ci>5^ v , B is the beam pattern, and iV^ rs crfd^ 
is the covariance matrix of the noise in the sky map. (Here Greek indices stand for pairs 
of indices (£m), with fx = t 2 + i + m + 1.) 
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The Karhunen-Loeve transform is a prescription for choosing the basis functions fa, or 
equivalently the elements of the matrix F. We choose the functions so that the likelihood 
in eq. (6) will have maximal rejection power for incorrect models. We therefore want 
the likelihood L to be, on average, as sharply peaked as possible about its maximum, 
or in other words, we want to choose Fij to maximize (-L"), where the primes denote 
derivatives with respect to some parameter in our theoretical model and the derivatives 
are evaluated at the maximum of L. This optimization problem reduces to a generalized 
eigenvalue problem: each row /, of F satisfies the equation 

Vaigf = Wnoisef- (8) 

Furthermore, the rows should be chosen to have the maximum eigenvalues A. Rows with 
small values of A probe mostly the distribution of the noise, with little sensitivity to the 
cosmological signal. They can therefore be omitted from the likelihood estimates with 
little loss of information. We have found that it is necessary to keep only the 400 most 
significant modes. 

Since we have no knowledge of the true monopole and dipole in the sky map, we 
marginalize over these modes. The peak value, width and location of this final marginalized 
likelihood, as a function of D 1: D[ and D", are the output of the fitting procedure. 

Now to find the normalization of any theory, one calculates the large-angle multipole 
moments and finds the quadratic which best describes their shape. Over the range —0.5 < 
D[ < 0.5 and —0.5 < D'{ < 3.5, the best-fitting amplitude and likelihood are given by the 
following analytic forms: 

10 n Ci = 0.8073 + 0.0395D; - 0.0193D" (9a) 

In L = 0.00697 + 1.523-Dj - 0.403D; 2 - 0.490D" - 0.0391D[D'{ + 0.00855D'/ 2 (96) 

The fitting formula for C\q has a worst-case error of 2% and an average error of 0.4% over 
this range; the corresponding numbers for L (not InL) are 7% and 1.7%. The uncertainty 
in Cio is approximately 15% for all models. Note that the COBE data prefer models 
with positive D[ and negative D". The likelihood reaches its maximum at the point 
(D'±, D'l) = (0.0,-3.0), which is beyond the range covered by the fitting formula. The 
peak likelihood for this model is 3.7 times the likelihood of a flat Harrison-Zel'dovich 
model. Fig. 2 shows L function of D[ and D'{. 

In Tables 1 and 2 we show the best fitting shape parameters for some flat, low-fio 
variants of the CDM model. The fit of a quadratic to these theories gives an error at 
the worst fit multipole (in the range £ = 2 to 30) of about 5%, with a typical error of 
< 2%, showing that such theories are well fit by quadratics over the range of scales probed 
by COBE. To quantify the error introduced by approximating the power spectrum as a 
quadratic, we computed likelihood curves for the worst-fitting model in our sample using 
both the true power spectrum and the quadratic approximation. The two curves differ by 
11% in peak likelihood and by 0.5% in normalization. 

Critical (fi = 1) CDM models with late reionization are also well fit by quadratic 
D(x). The most plausible, though not the only, ionization history in hierarchical models 
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of structure formation is standard recombination, followed by full ionization from some 
redshift z* until the present. The fully ionized phase is due (perhaps) to radiation from 
massive stars on scales which go non-linear early. (See Liddle & Lyth 1995 for further 
discussion.) We find that, over the range of t probed by COBE, models which have z* < 100 
are almost indistinguishable from models with no reionization, assuming standard big-bang 
nucleosynthesis values for Og. There is of course damping on degree scales (£ ~ 100), but 
little change in the spectrum at smaller I. Further, the relative normalization of the matter 
and radiation power spectra is the same as in models with standard recombination. For a 
CDM model with h = 0.5 and n = 1, the quadratic parameters are well fit by the formulae 



over the range 30 < z* < 110. 

The case of open CDM models is more complicated, since the Cg exhibit features at 
several scales. To add to the difficulty, there appear to be several different primordial spec- 
tra one can consider in open universe models. Some models based on inflationary phases 
(Lyth & Stewart 1990, Ratra & Peebles 1994, Bucher, et al. 1995) predict power spectra 
which show an increase in power near the curvature radius. All of these calculations make 
use of basis functions in which there is exponential damping of power above the curvature 
radius; however, this assumption can be relaxed (Lyth & Woszczyna 1995, Yamamoto 
et al. 1995). For further discussion of these issues see the appendix. The open models 
of Ratra & Peebles (1994) have already been fitted to the 2-year COBE data (Gorski et 
al. 1995a) and we will not duplicate the results here. We mention, however, that the Cg 
for such models can be fitted by cubics in log 10 1 to the same accuracy that the ACDM 
models can be fitted by quadratics. This increases the dimension of the parameter space 
and makes tabulating the results more difficult. We defer consideration of cubic fits until 
the situation with regard to open models is more settled. 

Once the 4-year COBE data becomes available, we hope a fitting formula similar to 
Eq. 9 (but which goes to sufficiently high order to encompass almost all theories), could be 
produced for the benefit of the astrophysics community. Such a fit, coded into a subroutine, 
would allow any theory to be quickly and accurately fitted to the COBE data. At present 
our simple quadratic fit is sufficient for a wide range of theories of current interest. 

The goodness of fit 

Statistical methods for using a data set like COBE to place constraints on models 
generally come in two varieties, Bayesian and frequentist. Most CMB work, including 
analyses of the COBE data as well as other experiments, has taken a Bayesian point of 
view. In the Bayesian approach, the probability of observing the actual data is computed 
for each model. This may be denoted p(D\M), meaning "the probability of the data given 
the model". One then assumes a "prior" probability distribution p(M) on the models 
and applies Bayes's theorem to produce a "posterior" probability distribution giving the 
likelihood p(M\D) of the various models given the data: 



D[ = 0.738 - 0.0307z* + 3.32 x 10" 4 ^ - 1.06 x 10" 6 ^ 
D'{ = 1.554 - 0.0483z* + 3.67 x 10" 4 ^ - 7.60 x 10 -7 ^ 



(10) 



p(M\D) oc p(D\M)p(M). 



(11) 
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The posterior probability distribution tells us how likely each member of our family of 
models is compared with any other member, which is what we would like to know. 

The Bayesian approach is perfectly adequate for assessing the relative merits of the 
various models under consideration: with this approach we can say, for example, that model 
A is 10 times more likely than model B. These models may only differ by a normalization 
or could be drawn from different cosmologies or structure formation scenarios. In some 
cases however we would like to assign absolute consistency probabilities to models. The 
Bayesian approach is not well suited to answering this sort of question, and the problem of 
assigning an absolute consistency probability to a model is best attacked with frequentist 
methods. 

In the frequentist approach, we choose some goodness-of-fit parameter rj, and compute 
its probability distribution over a hypothetical ensemble of realizations of the model. We 
then compute the value of rj corresponding to the real data, and determine the probability 
of finding a value of rj as extreme as the observed value in a random member of our 
ensemble. We take this probability to be a measure of the consistency of the data with 
the model: if the data does not occur often in realizations of the model, we say the model 
is "unlikely" given the data. 

Two points about this technique deserve emphasis. First, the consistency probabilities 
derived in this manner are conceptually quite distinct from Bayesian likelihoods. Bayesians 
and frequentists ask different questions of their data, and will therefore sometimes get 
different answers. We do expect that models which have low Bayesian likelihoods will 
in general have poor frequentist consistency probabilities; however, there is no generally 
applicable quantitative relation between the two. Second, it is clear that the success of the 
frequentist approach depends on choosing an appropriate goodness-of-fit parameter rj. For 
some classes of problems a standard choice is available; for the problem we consider below, 
we are unaware of such a standard. This is because the measurement of CMB fluctuations 
involves detecting extra noise on the sky. Thus the correlation function, or errors on the 
temperatures, which are used in the fit depend on the theory being tested, unlike normally 
examined cases of model fitting. 

We wish to assign frequentist consistency probabilities to various power spectra. The 
first goodness-of-fit parameter one might think of for this purpose is a simple x 2 , 



where x\ is the amplitude of the ith element in our eigenmode expansion, af is the variance 
predicted for Xi by our model, and M is the number of modes in the expansion. (In 
order to remove all sensitivity to the monopole and dipole, the eigenmodes /j should be 
orthogonalized with respect to these modes before the Xi are computed.) This parameter 
would be a natural choice if we wished to constrain the normalization of a model; however, 
our primary interest is in constraining the shape of the power spectrum, and this goodness- 
of-fit parameter is not well suited for this purpose. In fact, given any power spectrum 
Cg, we can choose a normalization that gives a x 2 that lies exactly at the median of its 




(12) 
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probability distribution, since o~i scales with the normalization of the theory. We would 
therefore conclude that for some normalization this model is a perfectly good fit regardless 
of the shape of the Cg. 

To focus on the power spectrum, let us consider quantities quadratic in the amplitude 
of the eigenmodes. There is a complication due to the presence of the galaxy in the COBE 
maps, which breaks the rotational symmetry of the COBE sky. This makes it difficult to 
define a rotationally symmetric quantity (like Cg) by summing over azimuthal variable m. 
However, this is only a technical complication and we can still define a measure of power 
by binning the squares of the mode amplitudes in bins that probe particular angular scales. 
We expand each eigenmode fi in spherical harmonics and compute an "effective £" probed 
by that mode by performing a weighted average over t with weights given by the squares 
of the coefficients of the expansion. (The modes are generally quite narrow in £-space, so 
the results are not sensitive to the exact method of computing the effective £.) We then 
sort the modes in order of increasing effective I (decreasing angular scale). As it happens, 
the result of this procedure is almost identical to sorting the modes in order of decreasing 
signal-to-noise eigenvalue. We then compute the quantities 



for 1 < % < M/K. We should choose the bin size K to be large enough to reduce the 
intrinsic width of the distribution of Zi to a reasonable level, yet small enough that the 
mode amplitudes in each bin probe similar angular scales. We have adopted K = 10 as a 
compromise between these two considerations. 

If our model is correct, then each Zi will be approximately K. If the model is incorrect, 
then some z\ will be too low, and others will be too high. For example, if our model has 
too little large-scale power, then the variances xj/aj will be greater than 1 for small j, 
and the first few Zi will tend to be larger than K. We can quantify this observation by 
defining the goodness-of-fit parameter 



The Zi for a Harrison-Zel'dovich spectrum and our worst-fitting model (model 4 of 
Table 3) are shown in Fig. 3. (In making Fig. 3 we chose the coarser bin size K = 20 
rather than K = 10, to reduce scatter in the points.) Note that with our definition the 
Zi only loosely correspond to Cg and depend on the theory. When Zi is larger than its 
expected value, one can conclude that the data have more power than the theory on the 
corresponding angular scale; however, there is no direct proportionality between Zi and 
the corresponding Cg. Each z% can be regarded as an estimator of the power spectrum 
of the signal and noise combined. For small i, z% samples mostly signal, while the noise 
dominates for large i. The value of Zi for large i therefore changes very little as the model 
parameters are varied, as can be seen in Fig. 3. 




(13) 



M/K 




(14) 
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Since the mode amplitudes Xi are in general correlated, it is not possible to compute 
analytically the probability distribution of r\. We must therefore resort to Monte Carlo 
simulations. For each of several models, we created 1,000 random sky maps. We added 
noise to each pixel by choosing independent gaussian random numbers with zero mean 
and standard deviations corresponding to the noise levels in the real data. We computed 
the parameter r\ for each map. We chose to simulate six different models. The first four 
were chosen to span a range of values of the Bayesian likelihood: we simulated (1) a 
Harrison-Zel'dovich spectrum; the models with the (2) highest and (3) lowest Bayesian 
likelihoods from our grid of quadratic power spectra; and (4) a model with an even lower 
Bayesian likelihood L = O.OILhz- In addition, we chose two models from Table 1 which 
have identical cosmological parameters (O = 0.1, h = 0.75, n = 0.85), except that (5) one 
has only scalar perturbations and one (6) includes tensors in the ratio Cj/Cf = 7(1 — n). 

Our simulation procedure fails to mimic the real COBE data in at least two ways. 
First, the assumption that the noise in different pixels is independent is not strictly true 
(Lineweaver et al. 1994). However, the correlations are quite weak and have been shown 
to have negligible effects in analyses similar to ours (Tegmark & Bunn 1995). Second, we 
have not attempted to model the removal of systematic effects from the data. However, 
we expect this to have little effect on the final results, since the removal of systematic 
effects primarily affects the low-£ multipoles, while proper treatment of the noise is more 
important for the high values of I where noise dominates. 

It is clear from Table 3 that low Bayesian likelihoods tend to correspond to poor 
frequentist consistency probabilities, as expected. Furthermore, those models with likeli- 
hoods of order unity are reasonable fits to the data. This is a very reassuring fact: it was 
perfectly possible a priori that all the models we have been considering would prove to be 
intrinsically poor fits to the data. Comparison with the consistency probabilities for mod- 
els 5 and 6 and a look at Fig. 4 allows one to calibrate the sensitivity of the COBE data 
to spectral shape information. We expect this to improve with the 4-year data, especially 
at higher £. 

The matter power spectrum 

The best normalization and the goodness of fit of the temperature fluctuations for a 
range of models are given by Eq. 9. Using these results to normalize the matter power 
spectrum from the CMB can present some complications. In the simplest picture, in 
which large-angle CMB anisotropies come purely from potential fluctuations on the last 
scattering surface, the relative normalization of the CMB and matter power spectrum 
today is straightforward (White, Scott & Silk 1994, Bunn et al. 1995). In the conventional 
notation where the radiation power spectrum is given by 



T ( i + n f L ) r ( 2 + ^ir) 



(15) 




the matter power spectrum for an Oq = n = 1 CDM universe is 



-2.5 x 10 16 A(A;//iMpc- 1 )T 2 (A;) (/iMpc -1 ) 3 . 



(16) 
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with A = 3C 2 /(47r). In models such as CDM this relation works quite well, as long as 
matter-radiation equality is sufficiently early (h is not too low). Even for h = 0.3 the 
relation works at the 4% level and if h = 1 the error is < 1%. See Bunn et al. (1995) for 
further discussion. 

For models with Qq < 1 the normalization is not so straightforward. Naively one 
would think that, for fixed CMB fluctuations at z = 1, 000, one would have smaller matter 
fluctuations today. This is because in an open or a flat model with a cosmological con- 
stant (Qa = 1 — ^o) 5 density perturbations stop growing once either the universe becomes 
curvature or cosmological constant dominated (respectively). Curvature domination oc- 
curs quite early, and the growth of density fluctuations dp/p = 5 in an open universe is 
suppressed (relative to an fi = 1 universe) by a factor Oq 6 . In a flat A model, the cosmo- 
logical constant dominates only at late times and so the growth suppression is a weaker 
function of the matter content: 5 oc 0[j 23 . This suppression of growth in an f!o < 1 uni- 
verse has often been cited as "evidence" that fio must be large — otherwise fluctuations 
could not have grown enough to form the structures we observe today. 

In fact there are several other effects which come into play when normalizing the mat- 
ter power spectrum to the COBE data in a low-fi model. The first is that, though the 
growth in such models is suppressed by Oq (p ~ 0.6 for open and 0.23 for A models; for a 
more general formula see Carroll et al. 1992), the potential fluctuations are proportional 
to Qq. Hence the CMB fluctuations are even more suppressed than are the density fluctu- 
ations! So for a fixed COBE normalization the matter fluctuations today are larger in a 
low-fio universe, and the cosmological constant model clearly has the most enhancement 
since the fluctuation growth is the least suppressed. In terms of the power spectrum, P(k), 
we expect for fixed COBE normalization that P(k) oc 5 2 oc _ , as has been pointed 
out by Efstathiou, Bond & White (1992). 

This potential suppression is not the only effect which occurs in Iow-Oq universes, 
although it is the largest. Due to the fact that the fluctuations stop growing (or in other 
words the potentials decay) at some epoch, there is another contribution to the large-angle 
CMB anisotropy measured by COBE. In addition to the redshift experienced while climbing 
out of potential wells on the last scattering surface, photons experience a cumulative energy 
change due to the decaying potentials as they travel to the observer. If the potentials are 
decaying, the blueshift of a photon falling into a potential well is not entirely canceled 
by a redshift when it climbs out. This leads to a net energy change, which accumulates 
along the photon path. This is often called the Integrated Sachs- Wolfe (ISW) effect, to 
distinguish it from the more commonly considered redshifting which has become known as 
the Sachs- Wolfe effect (both effects were considered in the paper of Sachs & Wolfe 1967). 
This ISW effect will operate most strongly on scales where the change of the potential is 
large over a wavelength, so preferentially on large angles (Kofman & Starobinsky 1985). 

In A models the ISW effect can change the relative normalization of the matter and 
radiation fluctuations at the 25% level for Q ~ 0.3 (see below). We show in Fig. 5 how 
these various effects on the inferred matter power spectrum normalization scale with Oq 
in a cosmological constant universe (the simplest case). We see the total power is slightly 
changed, for fixed C\q, because the shape of the Cg depend on Qq. This affects the goodness 
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of fit with the COBE data (see Bunn & Sugiyama 1995 and our Eq. 9b). The ratio of the 
large-scale matter normalization to C\q is changed by the ISW contribution to Cio, the 
change in the potentials and the growth of fluctuations from z = 1, 000 to the present. 
Over the range Oq = 0.1 to 0.5 one finds for an n = 1 spectrum with Cio = 10 -11 



almost independent of h. This can be compared with the scaling presented above. Also 
the epoch of matter-radiation equality is shifted, which changes the normalization on 
smaller scales for fixed large-scale P(k). Putting these effects together we show the RMS 
fluctuation on a scale 0.028 /iMpc -1 (see below) as a function of fi m Fig- 6. The sharp 
downturn at low VLq is due to a combination of the larger scale of matter-radiation equality, 
moving the break in the power spectrum to smaller k, and the photon drag on the baryons 
having an increased effect on fluctuation growth for large O^/Oo- For fio — 0.3 the shift 
in matter-radiation equality and the scaling of Eq. (17) roughly cancel, making A 2 much 
less sensitive to Qo than the individual contributions would suggest. 

We note here that the shape of the Cg for the tensor (gravitational wave) modes is 
largely independent of fi (Turner, White & Lidsey 1993). For this reason the radiation 
power spectrum of a A model with some tilt and a component of tensors can exhibit less 
curvature at I ~ 10 than a purely scalar power spectrum (see Fig. 4). Since in some 
inflationary models we expect a non-negligible tensor component (Davis et al. 1992, but 
see Liddle & Lyth 1992, Kolb & Vadas 1993) we have computed the tensor Cg following 
Crittenden et al. (1993) and give results both including and excluding a significant tensor 
contribution. Our results update those of Kofman, Gnedin & Bahcall (1993) who also 
considered tilted, ACDM models with a component of gravity waves. 

In open models, where curvature domination occurs early, much of the large-angle 
anisotropy comes from the ISW effect (Hu & Sugiyama 1995) so the matter-to-radiation 
normalization is even more complicated. For open models the dependence of P(k — > 0)/Cio 
on Q is not well fit by a power law, since the shape of the Cg on all scales depends on fi - 
Fig. 6 shows the normalization on smaller scales for a model with P(q) given by Eq. A3, 
where q 2 = k 2 /(—K) — 1 and K = H 2 (Qo — 1). [The Cg in this model will be similar to 
those in the inflationary models of Bucher et al. (1995) and Yamamoto et al. (1995). For 
a discussion of P(q) in open models see the appendix.] Notice that, for fixed h and CMB 
normalization, the open models predict a smaller amplitude for the matter fluctuations 
today than the A models. We can understand this as a consequence of the earlier onset 
of curvature domination than A domination and the consequently stronger suppression of 
fluctuation growth in open models. 

In Tables 1 and 2, and in Fig. 6, we show the normalization of the matter power 
spectrum for a range of models where the CMB normalization (C w ) is held fixed. We 
quote both the value of the RMS density fluctuation at 0.028/iMpc _1 (large-scale) and <t 8 
(small-scale). For comparison, Peacock & Dodds (1994) give 




= 1.14 x 10 6 O, 



-1.35 




(/i^Mpc) 4 



(17) 



A 2 (k = 0.028/iMpc" 1 ) 



(0.0087 ±0.0023) ( 



-0.3 




(18) 



dink 
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There are many determinations of a 8 ; we quote here those from Peacock & Dodds (1994) 



o-g = 0.75O, 



-0.15 

'0 



(19) 



and cluster abundances (White, Efstathiou & Frenk 1993) 



a 8 = 0.57O, 



-0.56 




(20) 



where the scaling with Qo in both cases refers to models with Qo+^a = 1- These values are 
consistent with those inferred from large-scale flows (Dekel 1994) and direct observations 
(Loveday et al. 1992). In Fig. 7 we compare these observations with the COBE-normalized 
values of cr 8 for a range of ACDM models. 

Conclusions 

The COBE data forms a unique and valuable resource for the study of inhomogeneous 
cosmology. To fully exploit this hard won information we need to go beyond methods of 
normalizing theories of structure formation which use only gross properties of the data 
(such as the RMS fluctuation). In this paper we have presented a model independent 
method of parameterizing the COBE data, and discussed the normalization of the radiation 
and matter power spectra for a range of theoretically interesting models. In addition we 
considered the question of the goodness of fit of some commonly adopted models to the 
data from two complementary statistical standpoints. 
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Appendix 



In this appendix we make some comments about the fluctuation spectrum in an open 
universe. The material is taken from the work of Lyth & Stewart (1990), Ratra & Pee- 
bles (1994, and references therein), Bucher et al. (1995) and Lyth & Woszczyna (1995). 
These papers are sophisticated and rigorous treatments of the subject, and consequently 
are somewhat lengthy and technical in parts. Here we try to give a flavor of the problem, 
building on the rigorous results of those works. We will proceed in historical order. 

We start by recalling a few points about inflationary cosmology. In an inflationary 
theory with fio = 1, quantum fluctuations during inflation give rise to density and poten- 
tial fluctuations today (see e.g. Olive 1990, Mukhanov et al. 1992, Kolb & Turner 1990, 
Linde 1990 and references therein). The amplitude of the fluctuations is set by the Hub- 
ble constant, H in ^ when the perturbation crosses out of the horizon, with larger scales 
(today) crossing the horizon earlier during inflation. If the potential of the inflaton (and 
hence l?inf) does not change very much during the time fluctuations on the relevant scales 
are produced (exponential inflation) one obtains a scale-invariant spectrum of fluctuations: 
5 2 oc k. In terms of potential fluctuations, which are the perturbations which enter the 
underlying metric and are in some sense more fundamental than the density perturbations, 
this corresponds to 5$ =constant oc H in f/mp\. (This spectrum is therefore described as 
scale invariant.) In Lyth & Stewart (1990) it was shown that, for scales which entered the 
horizon "early" when curvature could be neglected, 5<& =constant for Qq < 1 also. The 
transformation from this statement about the fluctuations in $ per logarithmic interval 
to one about P(k) proceeds in two steps. First, we relate the density perturbations to 
the potential fluctuations through Poisson's equation, which in an open universe reads 
(Mukhanov et al. 1992) 

(k 2 - 3K) 5$ k = Ana 2 p5 k (Al) 

where the curvature scale is K = Hq(Q — 1) as before. Second, we note that in an 
open universe, the eigenfunctions of the operator V 2 with eigenvalues k > y — K form 
a complete set (Lifshitz & Khalatnikov 1963, Harrison 1970, Abbott & Schaefer 1986). 
Thus we can expand all perturbations in term of these eigenfunctions. It is convenient to 
introduce the new variable q 2 = k 2 / (-K) — 1 which runs from to oo. [In order to obtain 
the most general gaussian random field in an open Universe, one must in general include 
the eigenfunctions with < k < \f—K] however, for fluctuations generated by inflation 
only modes with k > \J — K are excited (Lyth & Woszczyna 1995)]. Recall it is dQ^dk/k 
or P(k)d 3 k which gives the physical fluctuations, so we need to compare the power per 
logarithmic interval in k to the volume element in the two coordinates 

dk d 3 k d 3 q , , . 

~k = ~B = {A2) 

Writing k 2 — 3K oc q 2 + 4, and using the Poisson equation to translate from potential 
fluctuations (squared, per In k) to density fluctuations (squared, per d 3 q) we have 

(q 2 + 4) 2 
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which is the result of Lyth & Stewart (1990). Ratra & Peebles (1994) performed a calcu- 
lation of fluctuations from a linear potential in which they showed that this result extends 
to all q, not just those which entered the horizon while O ~ 1. 

In a recent paper, Bucher et al. (1995) consider an explicit model for open universe 
inflation. [A similar calculation has been done by Yamamoto et al. (1995)]. In this model 
the inflaton first gets trapped in a false minimum of the potential for some time. During 
this time the universe inflates exponentially. The quantum fluctuations in the zero point 
energy are (power-law) suppressed by the existence of a mass gap (the inflaton has a 
mass, since the potential has curvature at the minimum V((p) = Vq + m 2 <p 2 /2 + •■•). The 
inflaton then tunnels through the barrier in the standard semi-classical way (c.f. nuclear 
decay) and nucleates a bubble of fi — universe. As the potential rolls slowly from its 
post-tunneling value to the minimum of the potential more fluctuations are generated. 
The upshot of this, after a strenuous calculation, is that the potential fluctuations on 
small scales are <5$ =constant, as before. On larger scales, corresponding to earlier times 
during inflation, the potential fluctuations are either enhanced or reduced, depending on 
the value of m in the potential. For the value of m considered in Bucher et al. (1995) one 
finds an enhancement by an extra factor of q~ x on very large scales. It is worthwhile to 
stress however that the CV are relatively insensitive to q < 1 within reasonable variations 
in P(q). 

The issue of fluctuations in open-universe inflation is not settled. Further theoretical 
work is required to determine whether inflation makes a unique prediction for the fluctua- 
tion spectrum in an open Universe, or whether we must rely on experiments to distinguish 
among the different possibilities. Should the power spectra rise as q~ x or q~ 2 on large 
scales, as currently predicted, then models with fi between 0.1 and 0.3 will be disfavored 
by the COBE data. 
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Table 1: The shape of the radiation power spectrum in a ACDM model with h = 0.75. 
Also shown is the matter power spectrum normalization with the radiation normalized to 
Cio = 10 -11 . For the tilted models we show the results with (right columns) and without 
(left columns) a gravity wave component with Cj/Cf = 7(1 — n). 



14 





n 


rv 


7-1// 
D l 


A2 
A 


0"8 


rv 


D l 


A2 


C8 


0.50 


0.85 


-0.001 


0.657 


0.018 


0.80 


-0.123 


1.039 


0.010 


0.60 




0.90 


0.130 


0.698 


0.020 


0.89 


0.028 


0.937 


0.013 


0.72 




0.95 


0.262 


0.769 


0.023 


1.00 


0.196 


0.877 


0.018 


0.89 




l.UU 


U.oyo 


U.oDD 


U.UzD 


1 1 o 
l.lz 


u.oyo 


U.oDD 


U.U/D 


1 in 

l.lz 


0.40 


0.85 


-0.018 


0.807 


0.018 


0.67 


-0.130 


1.099 


0.010 


0.50 




0.90 


0.111 


0.851 


0.021 


0.75 


0.016 


1.021 


0.013 


0.60 




0.95 


0.240 


0.924 


0.024 


0.84 


0.178 


0.992 


0.019 


0.74 




1.00 


0.372 


1.025 


0.027 


0.94 


0.372 


1.025 


0.027 


0.94 


0.30 


0.85 


-0.043 


1.104 


0.017 


0.51 


-0.140 


1.222 


0.009 


0.36 




0.90 


0.084 


1.154 


0.020 


0.56 


-0.000 


1.188 


0.012 


0.44 




0.95 


0.212 


1.227 


0.022 


0.63 


0.155 


1.212 


0.017 


0.55 




1.00 


0.340 


1.327 


0.026 


0.70 


0.340 


1.327 


0.026 


0.70 


0.20 


0.85 


-0.071 


1.668 


0.013 


0.32 


-0.147 


1.436 


0.006 


0.22 




0.90 


0.055 


1.719 


0.015 


0.35 


-0.015 


1.471 


0.008 


0.27 




0.95 


0.180 


1.788 


0.017 


0.39 


0.131 


1.594 


0.012 


0.34 




1.00 


0.307 


1.883 


0.019 


0.44 


0.307 


1.883 


0.019 


0.44 


0.10 


0.85 


-0.029 


2.910 


0.003 


0.13 


-0.112 


1.806 


0.001 


0.08 




0.90 


0.096 


2.946 


0.004 


0.14 


0.018 


1.978 


0.002 


0.10 




0.95 


0.222 


2.998 


0.004 


0.15 


0.163 


2.326 


0.003 


0.13 




1.00 


0.347 


3.072 


0.005 


0.17 


0.347 


3.072 


0.005 


0.17 



Table 2: As in Table 1 but with h = 0.50. 
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Model 


D[ 


D'i 
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Consistency 
Probability 


1 


0.4 


-1.6 


3.29 


78.0% 


2 








1.00 


86.3% 


3 


-0.5 


3.5 


0.091 


93.0% 


4 


-2 





0.014 


97.2% 


5 


-0.342 


2.455 


0.186 


92.4% 


6 


-0.275 


1.561 


0.310 


91.4% 



Table 3: Bayesian and frequentist measures of goodness of fit for six cosmological models 
(see text). The shape parameters D[ and D" are defined in the text. L denotes the 
Bayesian peak likelihood, normalized so that a pure Harrison-Zel'dovich Sachs- Wolfe model 
has L = 1. The consistency probability is the percentage of simulated data sets for which 
the goodness-of-fit parameter rj defined in Eq. 14 is less than the value found for the real 
data. 
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Figure Captions 



Fig. 1: The value of the RMS power, as measured by cr(7°) (solid) and cr(10°) (dashed), 
for the best fitting, tilted Sachs- Wolfe models, as a function of n. The fact that the RMS 
fluctuation of the best fitting model depends on n shows that there is more information in 
the COBE data than just the RMS power. 

Fig. 2: The likelihood, L, as a function of the power spectrum shape parameters D[ and 
D'[. The contours range from L = 0.5 to L = 3 in steps of 0.5, where L = 1 corresponds to 
a flat spectrum D[ = D'{ = 0. For values of (D[, D'{) for which the power spectrum goes 
negative over the range 2 < I < 30 (lower right corner) we have set the likelihood to zero. 

Fig. 3: The data Zi (see Eq. 13) normalized to a Harrison-Zel'dovich spectrum (triangles) 
and our worst fitting model, model 4 of Table 3 (squares). The z% shown were computed 
with a bin size K = 20. The solid and dashed lines show the expectation value of each 
Zi and approximate one-sigma deviations from it. (The Zi are only approximately x 2 
distributed and only approximately uncorrelated, so these estimated deviations are not 
precise.) The effective t values probed by the different modes are indicated at the top of 
the figure. These quantities were computed by expanding each eigenmode fj in spherical 
harmonics and computing the centroid of the distribution in £-space, as described in the 
text. Note that with our definition the Zi only loosely correspond to Ci and depend on the 
theory. At higher i the modes are sampling mostly noise and therefore change very little 
as the model parameters are changed. 

Fig. 4: The power spectrum of fluctuations for a ACDM model with and without the tensor 
contribution. The models shown have h = 0.75, fi = 0.1 and n = 0.85, 0.90, 0.95 (bottom 
to top). The solid lines are scalar only, while the dashed lines have C<f/Cf = 7(1 — n). 
Note the tensor+scalar models have less curvature than the scalar only models, which 
makes them a better fitted to the COBE data. 

Fig. 5: An example of how the normalization of the matter power spectrum depends on Qq 
in a ACDM model (with n = 1). All quantities are normalized to their values at Qq — 0.5. 
The solid line is the RMS temperature fluctuation with C\q fixed. The dotted line shows 
the ratio of the large-scale matter normalization (lim fc ^ P(k)/k) to C w . The dashed lines 
show the effect of the shift in matter-radiation equality on the small-scale normalization erg, 
holding the large-scale normalization limfc^o P(k)/k fixed. We show two models: h = 0.75 
(upper) and h = 0.50 (lower) with fl B h 2 = 0.0125. 
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Fig. 6: The normalization, A 2 (0.028/iMpc _1 ), as a function of fi fc> r open CDM (dashed) 
and ACDM (solid) models normalized to C\q = 10 -11 . In both cases the upper curves are 
for h = 0.75 and the lower curves are for h = 0.50, both with Og/s 2 = 0.0125. 

Fig. 7: The small scale normalization as vs Qq in ACDM models, with normalization set 
by Eq. 9. The dashed lines assume all the contribution to the temperature fluctuations 
measured by COBE come from scalar perturbations; the dotted lines are scalars + tensors 
with Cf/Cf = 7(1 — n). Slope (n) increases from 0.85 to 1.00 in steps of 0.05 with lowest 
n being lowest cr 8 . The two solid lines are two observational determinations of a 8 , the top 
line from cluster abundances and the bottom line from large scale structure (i.e. Eqs. 19, 
20). 
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Erratum 



The COBE Normalization of CMB Anisotropies 
Martin White & Emory Bunn 
Ap. J. 450, 477 (1995) 

In tables 1 & 2 the quoted values of A 2 and as assume a tensor-to-scalar ratio 
/Cf = 7(1 — n). This is the correct lowest order (in an expansion in 1 — n) ex- 
pression for the tensor-to-scalar ratio for Qq = 1. However the projection from the /c-space 
inflationary prediction onto the quadrupole, i.e. 1 = 2 mode, has a dependence on Q\ and 
n, which was neglected in calculating the entries in the tables. The "correction factor", 
/ T ^ S (0 A , n), defined through 

Cf/Cf = 7 (1 - n) f T/s (n A ,n) 

can be well fit by 

/ t/s (Oa, n) = 0.97 - 0.58(1 - n) + 0.25fi A - [l - 1-1(1 -n) + 0.28(1 - n) 2 ] Vt\ . 

This expression includes both the full n dependence (i.e. beyond leading order in 1 — n) of 
power-law inflation, for which exact expressions are available, and the ISW contribution 
to Cf • To good approximation the values of A 2 quoted in the tables should be multiplied 

by 

l + 7(l-n) 
1 + 7(1 - n)f T/s 

and the values of cr 8 by the square root of this. For reasonable values of Q\ and n this 
amounts to a ^ 10% correction. 
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